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Abstract 

In supersymmetric theories, sneutrino-anti-sneutrino mixing can occur with 
the oscillation time ~ 0.01 ps corresponding the atmospheric neutrino mass 
scale ~ 0.05 eV. We explore the possibility of observing sneutrino oscillation 
phenomena and CP violation when R-parity violation explains the observed 
neutrino masses and mixing. It is shown for some parameter region in the 
bilinear model of R-parity violation that the tiny sneutrino mass splitting and 
time-dependent CP violating asymmetries could be measured in the future 
experiments if the tau sneutrino is the lightest supersymmetric particle. 



The atmospheric neutrino data from the Super-Kamiokande experiment strongly 
suggest the presence of sub-eV neutrino masses. As is commonly accepted, such tiny 
masses are most likely to be of the Majorana type breaking lepton number by two units 
(AL = 2). In supersymmetric theories where each neutrino species v is accompanied by 
a complex scalar v, similar lepton number violation should occur in the scalar (sneutrino) 
sector. The sneutrinos v and z>* generally have the mass terms, 

- £ = m\i>v* + i {m 2 M vv + /i.e.) , (1) 

where m| is the usual slepton mass of the order 100 GeV and m 2 M carries the same lepton 
number as the Majorana neutrino mass m v . As we will see, the soft 5-terms are the origin 
of the mass m 2 M in our case. This 'Majorana' sneutrino mass term generates mass splitting 
and mixing in the sneutrino- anti-sneutrino system [fj, |3|. One can easily find from Eq. (1) 
that the sneutrino mass-squared eigenvalues are m~ ± In general supersymmetric 

models, one finds m? M ~ m v my and thus the sneutrino mass splitting of the order of the 
neutrino mass; Amy ~ \m 2 M \/m v ~ m u . Similarly to the case of neutrinos, such a small 
mass difference could only be measured through the observation of sneutrino oscillation, 
providing a novel opportunity to probe AL = 2 lepton number violating phenomena in 
the future collider experiments ||. 

Mixing phenomena in the sneutrino-anti-sneutrino system are analogous to those in 
the B B system for which we refer the readers to the reviews in Ref. Extending 
the original investigation in Ref. ||, we wish to formulate CP violating effects in time- 
dependent sneutrino oscillation and apply them to a specific model of neutrino masses, 
namely, supersymmetric standard model with R-parity violation. For this, we assume the 
neutrino masses explaining the Atmospheric neutrino data with the largest neutrino mass 
m v . A ~ 0.05 eV and the large v^—Vt mixing []!]]. It is amusing to notice that the neutrino 
mass m U3 corresponds to the proper sneutrino oscillation time, t osc = Am^ 1 = ~ 
0.013 ps where L measures the distance between the sneutrino production and decay 
vertex and — is the boost factor of a few. Therefore, considering a spatial resolution of 
a few hundredths of ps ||, the atmospheric neutrino mass scale could be directly probed 
by observing time-dependent sneutrino oscillation in the future experiments. 
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For sneutrino oscillation to occur, it is required that the sneutrino life-time is longer 
than its oscillation time; 

Tc, < Amc, . 

In the supersymmetric see-saw model considered in Ref. the sneutrino mass splitting 
is given by Am^/m,, « 2(A + fi cot (3) /my where A is the trilinear soft mass and tan/? 
is the ratio of the vacuum expectation values of two Higgs bosons. With m v ~ 0.05 eV, 
one generically expects Am u < 1 eV. Then, the requirement Yy < 1 eV can be arranged 
in a limited region of parameter space under the assumption that the sneutrino has only 
three-body decay channels and the stau tr is the ordinary lightest supersymmetric particle 
(LSP) H. In this scheme, CP violating effects could arise if a nontrivial phase exists in 
the 'Majorana' sneutrino mass. 

Another interesting scheme for generating neutrino masses and mixing is to allow 
R-parity violation in the supersymmetric standard model J71. In this case, the LSP is 
unstable and decays through R-parity violating couplings which generate both neutrino 
masses and sneutrino mass splittings. Therefore, it is expected that R-parity violation 
can naturally lead to the suppressed decay rate and observable sneutrino oscillation phe- 
nomena if a sneutrino is the LSP. Whether or not a sneutrino can be the LSP depends 
on the models of supersymmetry breaking. In the constrained framework of minimal su- 
pergravity || or gauge-mediated supersymmetry breaking U, the sneutrino LSP may be 
obtained in a limited parameter space where the D-terms can make m~ smaller than rri~ R . 
But, such a possibility appears to be ruled out JL0| by current experiments on the invisible 
Z width, providing the limit m u > 44.7 GeV H. Therefore, we assume less constrained 
(non-universal) soft masses so that slepton doublets are lighter than slepton singlets and 
thus sneutrinos lighter than charged sleptons. Most favorable framework for this is SU(5) 
grand unification theory where the common soft mass of 5 multiplets is smaller than that 
of 10 multiplets or gauginos. Let us remark that our discussion can be complimentary to 
the case of a neutralino LSP in testing the R-parity violation model of neutrino masses 
and mixing [IT]. 

We begin our main discussion by considering the general sneutrino-fermion-fermion 
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couplings as follows: 
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(2) 



where Pl,r = (1 =F 7s)/2 and the index i refers to the lepton flavor. Neglecting mixing 
between different sneutrino flavors, each sneutrino species &i has the mass term (|T|) with 
complex m 2 M . The mass eigenstates denoted by two real fields, Dn and z> i2 , have the 
following decay widths: 

,2 
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Y, [K^K^e**" + c.c. 
/,/' 



(3) 
(4) 



where 4>m = Arg(m^). The time-evolution of the state identified as the sneutrino i>i or 
the anti-sneutrino v* at an initial time t = is given respectively by 



g+ (t)\%)+g4t)e«*»\fy 



(5) 



where g±(t) = \ exp(— ^T^t — imo a t)[l ± exp(|Ar^ + iAm^t)]. Then, the time 
dependent CP asymmetry for the sneutrino decay to the final state //' is 

nm^fw-Tww^ffi) 



Ahp(ff';t) = 



r(*i(*)->//0 + r(fi?(t) ->//') 

cos(Am Vi t)X{ f - 2 sin(Amp.t) lm(Y/ f ) 



(6) 



c(t)-s(t)Re(Y/ f ') 
where c(t) = {e AT ^ 2 + e~ Ar ^ 2 )/2, s(t) = {e AT ^ 2 - e - Ar M/ 2 ) and 



if 



\K{S\ 2 +\Kg\ 2 -(f^f>) 



Y; 



ff 



\K[ f i\ 2 +\K f R \{\ 2 + (f <-> /' 

K tr K ff + (/ <_> r) 



\K f 4\ 2 +\K f R j\ 2 + (f~f') 
For the 'flavor-specific' final state //' arising from the decay of Pj, that is, V{ — > ff' i>, 
we have = 1 and Y/^ = 0. Thus, the CP asymmetry in Eq. (||) becomes 

^ CP {ff ] t)=cos{Am,i) for f^f (I 
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when Ar,> « 0. This time-dependence could be measured to determine the sneutrino 
mass splitting as far clS Xi IS not too small. On the other hand, the time-integrated 
mixing probability ||, 

Xl ~ 2{x* + 1) W 
where = 2Am^ i /T i ) i and y,i = AT^./T^. can be used to determine the sneutrino mass 
splitting when Xj ~ 1. Recall that Xi can De measured by counting the 'same-sign' and 
'opposite-sign' lepton events, vp* — > lll'V and ViV* — > WW, analogous to the B system. 
As we will see, our model generically has j/j < 1. 

For the final state // which is shared by the decays of z/j and u*, we get X{* = and 
thus the time-dependent CP asymmetry (for yi w 0) becomes 

Acptfh t) = : §P , 2 sin(0 D - M sin Am,i (10) 
1 + \Pi\ 

where = \pi\e l< ^ D = K^* j K^ { . This is analogous to the CP asymmetry from I? and 
-B decays to CP eigenstates M. 

Let us now examine how the quantities discussed above arise in the R-parity violating 
supersymmetric standard model. For the sake of simplicity, we introduce only bilinear 
R-parity breaking terms in the superpotential W and the soft scalar potential V so f t of the 
super symmetry standard model; 

W 3 -toLiH-i, V soft 9 m 2 iH LiHl + B^Hi + h.c. (11) 

where the same notations Lj, Hi 2 are used for the lepton and Higgs superfields and their 
scalar components. This type of model is known to generate viable neutrino mass matrices 
explaining both the solar and the atmospheric neutrino data As is well-known, the 



bilinear terms in Eq. (|TT|) give rise to nonzero sneutrino vacuum expectation values 0; 
<H = (v*)/(Hl) = {mj H + /i/ij + B{t p) /ml. where tp = tan/? = (H®) / (H°) and /x is 
the supersymmetric Higgs mass parameter. Then the quantities /ij — //a, determine the 
tree-level neutrino mass matrix; 



m\ 3 = -M 2 Z ^4/F N (12) 



4 



where & = a* - m/ n and F N = M 1 M 2 /{M l c^ v + M 2 s 2 v ) - Mfsin2/3//z. Here M x and 
M2 denote the U(l) and SU(2) gaugino masses, respectively. Remember that only one 
neutrino, z/3, becomes massive from Eq. (|T2|), while the other two will get smaller masses 



from one-loop corrections. The mass matrix (|12|) fixes two neutrino mixing angles, #23 and 



#13, corresponding to the atmospheric neutrino and the reactor neutrino mixing angles, 
respectively, as follows flTTf: 



sin 2 2# 2 3 = 4|£ M | 2 |£ T | 2 , sin 2 2£ 13 = 4|£ e | 2 (l - |£ e | 2 ) (13) 

where £j = and |£| 2 = J2i |£j| 2 - Current experiments require ~ |£ T | for the large 

atmospheric neutrino mixing [[J and |£ e |/|£r| < 0.3 for the suppressed reactor neutrino 



oscillation v e —> v jliT [ 13] 



R-parity violation induces also nontrivial mixing between sneutrinos and neutral Higgs 
bosons. Furthermore, their CP-even and CP-odd parts mix together if general complex 
couplings are allowed. Therefore, we have to deal with a mass matrix of 10 neutral boson 
fields including the Goldstone mode. To do this, it is convenient to define the 'proper' 
sneutrino fields getting rid of the Goldstone mode by performing the see-saw rotation 
with the small angle dj. In this basis, one finds that the sneutrino-Higgs mixing term 
for each sneutrino generation is proportional to Bi — Ba>i [nj] . The next step is to rotate 
away these mixing terms to find the sneutrino mass splitting; 

Am 5j = 2m^M 2 z m\\r ]i \ 2 clsyF s (14) 

where rji = cii — Bi/B and F$ = {m 2 ,. — m 2 )(m?. — m 2 H ){m'~,. — m 2 ^). Here and m# 
denote the light and heavy CP-even Higgs boson masses, respectively and is the CP- 
odd Higgs boson mass which are defined in the R-parity conserving limit. Eq. ([14]) is 
consistent with the result of Ref. || for the CP-conserving case. In deriving the above 
result, we neglected the sneutrino flavor mixing induced by the neutrino flavor mixing. 
This is a good approximation as far as the mass differences between two sneutrino flavors 
are much larger than m„, — ^> m„, which is usually the case. Remark that the 
universality in B parameters, Bi/B = fii/fi, implies £j = rji which is not assumed in this 
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paper. Combining Eqs. fll2|) and (0), one finds 

~ ^/m^XMVlei 2 ) (15) 

for <C rriA- Therefore, we get Amp i ~ m Uz as expected for Fn ~ m,>. and j^l ~ |£|. 
Here we note that the sneutrino oscillation time l/Am^. can be made larger than l/m V3 ~ 
0.013 ps for |?7j| < |£|. From the discussion below Eq. (|T3|) implying |r/ e | ~ |£ e | <C |£|, one 
finds that the electron sneutrino v e will generically have the largest oscillation time. But, 
we will see that the u e decay to charged leptons have too small branching fraction to be 
observed. 

The sneutrino couplings to light fermions / and /' arise from the neutrino-neutralino, 
the charged lepton-chargino, and the slepton-Higgs mixing. In the bilinear model under 
consideration, one finds the coupling constants K L R defined in Eq. (g) as follows: 

KZ k = ^^tcp, (16) 



K L,i = m-Wt— + kT —\ , K Li = — — — 
M w Cp Cp M w Cp 

jf m T u* 2 9* 2 

K R,i = T7 [Oij — + Ojr — \ , K 



jf _ mr u* 2 9* 2 bb m b 9* 2 
M w Cp Cp M w Cp 



where the coefficients u and 6 are given by 



_ -fiM 2 & 2/x(// + M 2 tp)M^j k 2 

Ukl ~ »M 2 - M* w s 2 p + ° fc ' Uk2 ~ {/M - M* w s 2 pf ^ 
On = -at ~ s 2 p[ml t - m\m\ - (mj. + m 2 A c 2 p)M 2 z s 2 p\m 2 A 'q i / F s 

9 i2 = - C ypM 2 z [ml - mlc^mlrji/Fs . (17) 

In Eq. (|T6D, we neglected the terms proportional to lighter quark and lepton masses. In 
the presence of trilinear R-parity violating terms, the couplings (fIB|) will get additional 



contributions leaving Eqs. ( 12Q and ( |14|) unchanged. 



We are ready to discuss how the observable sneutrino mixing phenomena occur in 
our scheme. From Eqs. ( |TE{ ) and flT7|), one can see that the dominant decay channels are 
Vi — > vv and vi —> bb, rf for low and large tan (3, respectively, if the R-parity violating 
parameters 0Oi,0{, are of the same order. For the latter decay channels, one also finds that 
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\Kl\ is larger than \Kr\ for the parameters Mz < M 2 , p and large tan/3. This shows 
that we generally have y = AT/T <C 1 as the decay rate difference Ar gets non- vanishing 
contributions from Kjj K{/ where f = r, b. Since we are interested in the region where 
Yr, < Amr,, we can put c(t) = 1 and s(t) = in Eq. (g). Let us now make an estimate of 
the quantity x^. Taking T„. dominated by the couplings K v ^ i or KJ^, we get 





Vi 
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(M Z M W ) 2 
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(m^m T tp) 2 
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(18) 



which shows Xi 3> 1 for small tan/3. With = Mz = -Fat/2 and \r)i\ = |£| = \9a\, one 
gets Xi ~ 1 with tan /3 ~ 20 for which the time-integrated quantities can be well measured. 
On the other hand, for low tan /3, the time-dependent observables may be measured to 
determine Am^ or CP asymmetry. If tan/? is too large (xj 1), no oscillation effect 
occurs. The quantity p, determining the magnitude of CP asymmetry is given by 

Pi = i2 /9ii or (S iT u i2 + 9 i2 )/(5i T u il + 9 a ) (19) 

for the shared final states bb or rf, respectively. As can be inferred from previous discus- 
sions, the magnitude of pi is usually smaller than 1 and drops down with growing tan /3. 
This should be compared with the quantity Xi which shows the similar behavior as above. 
Thus, we typically get \pi\ <C 1 in the region where Xi ~ 1. For the measurement of the 
CP asymmetries in the shared final states, one should determine whether the decaying 
particle is a sneutrino or an anti-sneutrino. Obviously, this can be done by looking at the 
flavor-specific decays of the other particle, for instance, u a —> af </- v* a and v T —>■ ra *f- v* 
where a = e, p. Therefore, the branching fractions for the flavor-specific as well as the 
shared final states should be large enough. To get a reference value, let us consider a 
future e + e~ linear collider with an integrated luminosity 1000 fb _1 at y/s = 500 GeV. For 
the cross-section of the sneutrino-anti-sneutrino pair production ~ 50 fb with m 5 ~ 100 
GeV |p.5| , there will be about 5 x 10 4 events. Requiring the multiplied branching ratio 
to be larger than 1 %, one will get about 500 samples to study CP violation. Here, we 
note that such a value can be hardly obtained for the electron and muon sneutrinos. This 
is because the decay v a — > af (a = e, p) comes from the coupling K^ T a which is sup- 
pressed by the factor m T J tan/3. Considering the ratio of K^ a to K^ a and KJJ a , we get 
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the branching fraction, BR(af) = \K^ T a / K^ a \ 2 {m T tpM w / jj,M z ) 2 for low tan/?, and 
BR(af) = \K^ a / K T jJ a \ 2 « (M^ v /t^fiM 2 ) 2 for large tan/5, which shows that the branching 
ratio is usually smaller than 1 % for My/ < M 2 ~ Fn ~ /i. 

However, the situation can be different for the tau sneutrino which allows BR(r/x) ~ 
BR(rf) ~ BR(66) ^> BR(w) for |a; Ml | ~ |a; Tl | ~ |0 Tl | and relatively large tan/3, as 
the relevant coupling IC£ a T has no such suppression. To show explicitly the sneutrino 
oscillation and corresponding CP violation effects, we choose the following typical set of 
mass parameters reference: 

M 2 = n = 2m L , m A = 2.5m L , m h = 115 GeV (20) 

where m L denotes the slepton doublet soft mass. Taking m L = 100 (200) GeV, the 



sneutrino masses are, rxiy = Jm 2 L + Mf c 2j g/2 = 82, 77 and 76 (192, 189 and 189) GeV for 
tan /3 = 3, 15 and 30, respectively. Concerning the gaugino masses, the unification relation 
is assumed; M 2 = 2M\. In our bilinear model, there are three types of R-parity violating 
parameters /ij, Bi and m 2 H without assuming the universality. For our calculation, we 
will trade these parameters with rji and eij. Then, we put ^ = £ T — 3^ e , consistently 
with the experimental data as mentioned before and the value of £ T is normalized to yield 
m U3 = 0.05 eV in accordance with Eq. (0). In fact, our results are insensitive to the ratio 
£ e /£ T . Note that small neutrino masses requires very small R-parity violating numbers, 
£iC/3 ~ 10~ 6 , which could be a consequence of some flavor symmetry explaining quark and 
lepton Yukawa hierarchies [16[ . 



To find the favorable parameter space, we consider the following two regions: (a) 
\£,t\ ~ \Vt\ ~ | Or | an d (b) |£ r | ^> \r] T \ ~ \a T \. As discussed before, the latter region 
is chosen to give larger |£ r | and thus larger oscillation time. In the first region, the 
oscillation time becomes too small for the time-dependence to be observed. In Table A 
and B, we present the branching ratios and various oscillation parameters for the tau 
sneutrino LSP taking the following specific numbers for the two regions: 

(A) £ r = -r) T = -2a T , (B) £ T = 6 Vt = -6a T . (21) 

From Table A, one finds that x becomes of order one for tan/3 ~ 30 and the proper 



oscillation time t osc is very small. As the branching ratios for i> T — > rf , r/2 become 40- 
50% for large tan (3, one could look for the same-sign lepton signal to measure the time- 
integrated mixing probability x- 111 this case, the CP asymmetries p are smaller than 1 % 
and thus can be hardly measured. As alluded before, the set (B) gives observable time- 
dependence. As can be seen from Table B, p becomes smaller and the branching ratios 
larger as tan/5 grows. The multiplied branching ratios for the leptonic modes becomes 
of order 10% for intermediate to large tan/5 and thus the sneutrino mass splitting can 
be directly obtained by tracing time-dependent oscillation for the final state r/2 in this 
region. Furthermore, the time-dependent CP asymmetry for the final state rf can be a 
few percent assuming the maximal CP phase | sin(0£) — 0m) | ~ 1; which could be within 
the future experimental reach. In fact, the quantity p can be made larger for |£| ^> \r)i\ 
and lower tan/?. However, restricting ourselves to the region of |£| < 10|7/j|, we find that p 
can be maximally a few %. Therefore, a better luminosity than mentioned before will be 
needed to cover more parameter space of our model. We note here that one hardly gets 
large CP asymmetry for v T — > bb as the corresponding p is proportional to t^ 2 . In both 
cases of (A) and (B), one finds that smaller w,l is favored for the observation of various 
sneutrino oscillation observables. 

Let us finally comment on the case of the universal soft parameters giving & « r]i -C a^. 
From Eq. (|16|) and (]17[) , one can see that the sneutrino decays dominantly to bb and the 
branching fraction for rp is mj/3mj ~ 5% with much more suppressed rate for rf. Thus, 
almost no oscillation effects can be observed in this case. 

In conclusion, R-parity violation can lead not only to the realistic neutrino masses 
and mixing, but also to observable sneutrino-anti-sneutrino mixing phenomena. This can 
occur if the tau sneutrino is the LSP. Among various observable quantities in the sneutrino 
oscillation, a certain observable will only be within experimental reach given parameter 
region of the bilinear model of R-parity violation. For the less fine-tuned bilinear R-parity 
violating parameters, |£| ~ \r]\, only time-integrated quantities can be probed, whereas 
the time-dependence and CP asymmetries in the sneutrino oscillation can be observed 
for the region where |£| 3> \rj\. In the former region, the sneutrino mass splitting can be 
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determined by measuring the time-integrated mixing probability for large tan (3 where x 
becomes order one. In the latter region, the mass splitting can be determined through 
time-dependent oscillation even for lower tan f3. Finally, if the sneutrino mass is small 
(my < 100 GeV) and tan f3 is in the intermediate region, the CP asymmetry in the decay 
v T — > rf can reach the level of a few percent which could be measured in the future 
experiments. 
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m L /GeV 


tan j3 


vv 


TT 


T/J, 


bb 


X 


y 


p{rf) 


p{bb) 






3 


0.68 


0.03 


0.07 


0.22 


13035 


0.08 


0.037 


0.20 


0.008 


100 


15 


0.007 


0.38 


0.47 


0.09 


197 


0.02 


0.013 


0.04 


0.007 




30 


5xl0" 4 


0.40 


0.48 


0.08 


13 


0.007 


0.007 


0.012 


0.007 




3 


0.17 


0.21 


0.084 


0.53 


12670 


0.062 


0.031 


0.047 


0.0048 


200 


15 


4xl0" 4 


0.24 


0.092 


0.66 


35 


0.005 


0.003 


0.003 


0.0037 




30 


2xl0~ 5 


0.24 


0.092 


0.66 


2.2 


0.002 


0.001 


0.0007 


0.0037 



Table A: The branching ratios of the tau sneutrino decay into vi>, tt, t]x and bb are shown 
together with the sizes of the oscillation parameters x, y and p defined in Eqs. and 
(JT0|) , respectively. The choice of R-parity conserving parameters is described in Eq. ( |20| ) 
and the set (A) in Eq. (|2l|) is used for the R-parity violating parameters. 



m L /GeV 


tan p 


vv 


TT 


Tp 


bb 


X 
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P (tt) 


p(bb) 


tosc/ PS 




3 


0.87 


0.05 


0.06 


0.01 


468 


0.01 


0.16 


0.14 


0.30 


100 


15 


0.011 


0.27 


0.41 


0.28 


7.7 


0.011 


0.026 


0.005 


0.25 




30 


7xl0" 4 


0.27 


0.41 


0.29 


0.50 


0.006 


0.012 


0.001 


0.24 




3 


0.48 


0.06 


0.14 


0.30 


985 


0.015 


0.034 


0.017 


0.17 


200 


15 


2xl0~ 3 


0.11 


0.26 


0.60 


4.6 


0.003 


0.007 


0.001 


0.13 




30 


lxlO" 4 


0.11 


0.26 


0.60 


0.29 


0.001 


0.004 


0.0003 


0.13 



Table B: The same as above but with the set (B) in Eq. (fTTl). 
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